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We present direct numerical simulations of dynamo action in a forced Roberts flow. The behavior 
of the dynamo is followed as the mechanical Reynolds number is increased, starting from the laminar 
case until a turbulent regime is reached. The critical magnetic Reynolds for dynamo action is found, 
and in the turbulent flow it is observed to be nearly independent on the magnetic Prandtl number in 
the range from ~ 0.3 to ~ 0.1. Also the dependence of this threshold with the amount of mechanical 
helicity in the flow is studied. For the different regimes found, the configuration of the magnetic 
and velocity fields in the saturated steady state are discussed. 
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I. INTRODUCTION 

In a previous publication [J , a driven turbulent mag- 
netohydrodynamic (MHD) dynamo was studied numer- 
ically, within the framework of rectangular periodic 
boundary conditions. The emphasis was on the dynamo's 
behavior as the magnetic Prandtl number Pm (ratio of 
kinematic viscosity to magnetic diffusivity) was lowered. 
As Pm is lowered at fixed viscosity, the magnetofluid 
becomes more resistive than it is viscous, and it is in- 
tuitively apparent that magnetic fields will be harder 
to excite by mechanical motions. The principal result 
displayed in Ref. fIJ was a curve of critical magnetic 
Reynolds number, R C M , as a function of P M X , at fixed 
kinetic energy. The (turbulent) kinetic energy was the 
result of an external mechanical forcing of the Taylor- 
Green type (hereafter, "TG"), a geometry well known to 
be efficient at the rapid generation of small scales in the 
fluid flow 0. Ref. Q contains a lengthy bibliography of 
its antecedents, not all of which will be listed again here. 

The TG geometry injects no net mechanical helicity 
into the flow. In the long history of the dynamo problem, 
mechanical helicity has been seen often to be an impor- 
tant ingredient for dynamo action, and it is the intent of 
this present paper to consider a helically-forced dynamo 
in the same spirit as in Ref. fIJ , to see what changes oc- 
cur relative to the TG flow, further properties of which 
were displayed in a subsequent astrophysical paper 

A natural candidate for a highly helical velocity field 
is what has come to be called the "Roberts flow" 0, || . 
This flow shares some similarities with the dynamo ex- 
periments of Riga and Karlsruhe 0, 01 • In a pioneering 
paper || , Feudel et al. characterized mathematically var- 
ious magnetic-field-gencrating instabilities that a forced 
Roberts flow can experience. The present paper expands 
these investigations, while discussing numerical simula- 
tion results for magnetic excitations in the mechanically 
turbulent regime, with an emphasis on the nonlinearly- 
saturated magnetic field configuration. As in Ref. [s|, we 
will force the system at nearly the largest scale available 
in the periodic domain. As a result, magnetic fields will 
be only amplified at scales smaller than the energy con- 



taining scale of the flow. The behavior of the large-scale 
dynamo (i.e. when magnetic perturbations are amplified 
at scales larger than the energy containing eddies) as Pm 
is varied will be studied in a future work. 

Section [H] lays out the dynamical equations and defi- 
nitions and describes the methodology to be used in the 
numerical study. Section 1 1 i 1 1 presents results and com- 
pares some of them with the corresponding TG results. 
Section llVl summarizes and discusses what has been pre- 
sented, and points in directions that we believe the results 
suggest. 

II. MATHEMATICAL FRAMEWORK AND 
METHODOLOGY 

In a familiar set of dimensionless ( "Alfvenic" ) units the 
equations of magnetohydrodynamics to be solved are: 

— +v Vv = -W+j x B + j/V 2 v + f, (1) 
at 

— — h v ■ VB = B ■ Vv + 7/V 2 B, (2) 
at 

with V • v = 0, V • B = 0. v is the velocity field, 
regarded as incompressible (low Mach number). B is 
the magnetic field, related to the electric current den- 
sity jbyVxB=j. V is the normalized pressure-to- 
density ratio, obtained by solving the Poisson equation 
for it that results from taking the divergence of Eq. (JJJ 
and using the incomprcssibility condition V ■ v = 0. In 
these units, the viscosity v and magnetic diffusivity r\ can 
be regarded as the reciprocals of mechanical Reynolds 
numbers and magnetic Reynolds numbers respectively, 
where these dimensionless numbers in laboratory units 
are Ry = LU/vi a b, Rv — LU/viab- Here U is a typical 
turbulent flow speed (the r.m.s. velocity in the following 
sections), L is a length scale associated with its spatial 
variation (the integral length scale of the flow), and uiab, 
Viab are kinematic viscosity and magnetic diffusivity, re- 
spectively, expressed in dimensional units. The external 
forcing function f is to be chosen to supply kinetic energy 
and kinetic helicity and to maintain the velocity field v. 
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For f , we choose in this case the Roberts flow 0, Q : 
f = -vV 2 v R = 2vw R (3) 

where 

V -R = (.g sin x cos y,—g cos xsiny, 2/ sin x sin y) . (4) 

The coefficients / and g are arbitrary and their ratio 
determines the extent to which the flow excited will be 
helical. The ratio / = g/V2 is maximally helical for a 
given kinetic energy, and the case f/g — *■ is a (two- 
dimensional) non-helical excitation. We have concen- 
trated primarily upon the cases / = g (following Feudcl 
et al. 0) and / = g/y/2. No dyn amo can be expected 
unless \f/g\ > 0. 

We impose rectangular periodic boundary conditions 
throughout, using a three-dimensional periodic box of 
edge 2ir, so that the fundamental wavenumber has mag- 
nitude 1. All fields arc expanded as Fourier series, such 
as 

v = v(x, t) = ^ v(k, t) exp(ak ■ x) (5) 

k 

with k • v(k, t) = 0. The Fourier series are truncated at a 
maximum wavenumber k max that is adequate to resolve 
the smallest scales in the spectra. The method used is the 
by-now familiar Orzag-Patterson pseudospectral method 
0- EJ [HI ■ The details of the parallel implementations of 
the Fast Fourier Transform can be found in Ref . 01 ■ 

The forcing function Q injects mechanical energy at 
a wavenumber |k| = y/2, which leaves very little room 
in the spectrum for any back-transfer of helicity (|k| = 1 
is the only possibility). The phenomena observed will 
therefore be well-separated from those where an "inverse 
cascade" of magnetic helicity is expected to be involved. 
Rather, a question that can be answered (in the affirma- 
tive, it will turn out) is, To what extent does the presence 
of mechanical helicity in the large scales make it easier to 
excite magnetic fields through turbulent dynamo action? 

Equations © and |0J define a steady state solution 
of Eqs. © and @, with B = 0. It is to be expected 
that for large enough v and 77, this solution will be stable. 
As the transport coefficients are decreased, it will be the 
case that the flow of Eq. (0} can become unstable, cither 
purely mechanically as an unstable Navier-Stokcs flow, 
or magnetically as a dynamo, or as some combination of 
these. Thus rather complex scenarios can be imagined as 
either of the Reynolds numbers is raised. 

In the following, the emphasis will be upon discov- 
ering thresholds in Rm at which dynamo behavior will 
set in as Ry is raised, then further computing the non- 
linear regime and saturation of the magnetic excitations 
once it does. The "growth rate" a can be defined as 
a = dhx{E M )/dt, where E M = £ k l B ( k > *)IV 2 is the 
total magnetic energy. The appearance of a positive a 
for initially very small Em is taken to define the critical 
magnetic Reynolds number R C M for the onset of dynamo 
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FIG. 1: (Color online) Critical magnetic Reynolds R C M as a 
function of Rv for different Roberts flows (thick lines) : f = g 
(<>)> / = g/V% (□), / = 3/0.77 (A). The dark (red) area cor- 
responds to the region where the Roberts flow is hydrodynam- 
ically stable. For a comparison with Ref. 0, our Reynolds 
numbers should be divided by 2n. The light (orange) area 
corresponds to the region of hydrodynamic oscillations, while 
the white area corresponds to the turbulent regime. The thin 
lines connected with crosses are shown as a reference and 
correspond to the threshold for dynamo instability in Taylor- 
Green flow 0: DNS (solid line) and a-model (dashed line). 



action, a is typically expressed in units of the reciprocal 
of the large-scale eddy turnover time L/U where U is the 

1/2 

r.m.s. velocity (U = (u 2 ) , and the brackets denote 
spatial average), and L is the integral length scale, 

L = 2^^fc- 1 |u(k,t)| 2 /^>(k,i)| 2 . (6) 

k ' k 

In the next Section, we describe the results of the 
computations for both the "kinematic dynamo" regime 
[where j x B is negligible in Eq. (J)] , and for full MHD 
where the Lorentz force modifies the flow. 



III. DYNAMO REGIMES FOR THE ROBERTS 
FLOW 

We introduce the results for the Roberts flow through a 
plot of the threshold values of critical magnetic Reynolds 
number R C M vs. mechanical Reynolds number Ry (Fig. 
^1 . All Reynolds numbers have been computed using the 
integral scale for the velocity field [Eq. averaged 
over the duration of the steady state in the hydrodynamic 
simulation. For this same time interval, an overall nor- 
malization factor has been multiplied by Eq. J3} to make 
the r.m.s. velocity U turn out to have a time-averaged 
value of about 1. 

Figure n contains considerable information. There are 
basically three qualitative behaviors exhibited for differ- 
ent Ry, indicated by the (colored) background shading. 



3 



For Ry < 100, the laminar Roberts flow is hydrodynami- 
cally steady-state and laminar, but dynamo action is still 
possible for large enough Rm- For 100 < Ry < 1000, 
Roberts flow treated purely hydrodynamically is tempo- 
rally periodic but not turbulent. For Ry > 1000, the 
Roberts flow develops a full turbulent spectrum hydrody- 
namically. In all three regimes, dynamo action is exhib- 
ited, but is different in the three regimes. The laminar 
regime was extensively studied in Ref. Q. Our defi- 
nitions for the Reynolds numbers are different, but the 
results displayed in Figs. ^and[21are consistent with pre- 
vious results in the range Pm = [0.5, 1] if our Reynolds 
numbers are divided by 2ir (corresponding approximately 
to the integral scale of the laminar flow). 

The threshold curve connecting diamonds (O) is for 
the Roberts flow with / = g (helical, but not maximally 
so). The segment connecting squares (□) is for / = g/v2 
(maximal helicity). The segment connected by triangles 
(A) is for / = g/0.77, a less helical flow than / = g. The 
threshold curve connecting crosses (x) is the threshold 
curve for the Taylor-Green (TG) flow from Ref. 0|. All 
of these are direct numerical simulation (DNS) results. 
(We regard the fact that the Taylor-Green curve and the 
Roberts flow curve with / = g have a common region 
above Ry ~ 100 to be coincidental). The curve connect- 
ing crosses ( x ) with a dashed line is the result from Ref. 
[ij for the "a-model" , or Lagrangian averaged model, of 
MHD. 

Noteworthy in Fig. ^ is the qualitative similarity of 
the behavior of the threshold curve between the Roberts 
flow and the TG results from Ref. a sharp rise in 
R C M with the increase in the degree of turbulence in the 
velocity field, followed by a plateau in which further in- 
creases in Ry show little effect upon R C M - It must be 
kept in mind that for both situations, the amplitude of 
the forcing field f is being adjusted so that U and the to- 
tal kinetic energy remain approximately constant, even 
though Ry is increasing. Whether such a procedure cor- 
responds to a physical driving force must be decided on 
a case-by-case basis. 

Figure|2shows the threshold curve for the Roberts flow 
with / = g as a function of the inverse of the magnetic 
Prandtl number, P M . This curve shares some similari- 
ties with the TG flow, but also important differences. As 
in Ref. Q], between the laminar and turbulent regimes 
a sharp increase in R C M is observed. Also, in the tur- 
bulent flow R C M seems to be independent of the value of 
the magnetic Prandtl number. But while the TG force 
is not a solution of the Euler equations and was designed 
to generate smaller and smaller scale fluctuations as the 
Reynolds number Ry is increased, the Roberts flow goes 
through several instabilities as Ry is varied. As a result, 
the threshold for dynamo action in the Rm vs. P^ 1 plane 
is double- valued. For a given value of P M X two values of 
R° M exist according to the hydrodynamic state of the hy- 
drodynamic system, (e.g. laminar, periodic, or turbulent 
flow). 

Figure |21 is a plot of the kinetic energy spectra for the 
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FIG. 2: Critical magnetic Reynolds R C M as a function of P^j 
for the Roberts flow with / = g (thick lines). The solid line 
corresponds to the laminar regime, the dashed line to the pe- 
riodic flow, and the dotted line to the turbulent regime. (The 
double- valuedness results from the effects of two different val- 
ues of Ry-) 
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FIG. 3: Kinetic energy spectra as a function of Rv- The 
Kolmogorov's spectrum is showed as a reference. 

values of Ry shown in Fig. ^ for / = g, normalized so 
that Ey(k = 1) is unity for all cases. This is done to 
display the gradual widening of the spectrum as Ry in- 
creases. Figurc0]shows corresponding magnetic spectra, 
normalized somewhat differently: the energy contained 
in the interval 1 < k < 5 is the same in all cases. This 
is done to emphasize the fact that the peak in the mag- 
netic energy spectrum migrates to higher values as Ry 
increases: the excited magnetic field develops more and 
more small-scale features. This may be related to the fact 
that because the forcing occurs at such low wavenumbcrs, 
inverse magnetic helicity cascades are effectively ruled 
out. 

Figure [5] shows how the thresholds (a — 0) for the R C M 
curves were calculated. For small initial Em, broadly 
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FIG. 4: Magnetic energy spectra during the kinematic regime, 
for different values of Rv- The values of Rm for each curve 
correspond to the smallest value for which dynamo action was 
observed (see Fig. 
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FIG. 5: Growth rates as a function of Rm- Each line corre- 
sponds to several simulations at constant Rv (fixed v), and 
each point in the line indicates the exponential growth (or 
decay) rate at a fixed value of Rm. The point where each 
curve crosses a = gives the threshold R C M for dynamo in- 
stability. R v = 62.7 (□), Rv = 125.5 (x), R v = 416.2 
(A), R v = 974.4 (O), R v = 1108.1 (*), R v = 1327.7 (+), 
Rv = 1912.5 (o), and R v = 3276.6 (v)- 



distributed over k, r\ was gradually decreased in steps 
to raise Rm in the same kinetic setting until a value of 
a ~ was identified. That provides a single point on 
such curves as those in Fig. ^ 

Each simulation at a fixed value of v and r\ (or fixed 
Rv and Rm) was extended for at least 100 large-scale 
turnover times to rule out turbulent fluctuations and ob- 
tain a good fit to the exponential growth. All the sim- 
ulations were well-resolved and satisfied the condition 
ku/kmax < 1, where fc„ = (e/^ 3 ) 1 / 4 is the Kolmogorov 
lengthscale, e is the energy injection rate, k max = N/3 



FIG. 6: (Color online) Time history of the total kinetic 
[thick (blue) lines] and magnetic energy (thin lines) in dy- 
namo simulations. The dashed lines correspond to Rv = 62.7 
and Rm — 78.4 (laminar flow), while the solid lines are for 
Rv = Rm = 416.2. The shaded region indicates the period of 
time when the flow is oscillating in this simulation. The inset 
shows the time history for a turbulent run with Rv — 3276.6 
and R M = 1092.2. 



is the largest resolved wavenumber, and N is the lin- 
ear resolution of the simulation. When this condition 
was not satisfied, the resolution N was increased, from 
N = 64 until reaching the maximum spatial resolution 
in this work of 256 grid points in each direction, and a 
maximum mechanical Reynolds of Ry = 3276.6. 

Figure[S]illustrates an interesting behavior that occurs 
when there is a transition from the laminar to the peri- 
odic regime of the Roberts flow (/ = g). Figure shows 
the evolution of total kinetic energies Ey and magnetic 
energies E M for Rv = 62.7 and R v = 416.2. The flat 
part of the kinetic [thick (blue)] curve for Ry = 416.2 is 
characterized by small periodic oscillations too small to 
see on the logarithmic plot (they will be shown in Fig. [7J. 
Meanwhile, the Em curve of magnetic energy is growing, 
somewhat irregularly. Rather suddenly, at about t = 70, 
Ey drops by more than a factor of 2 (see Fig. 0) , and by 
t « 300 the magnetic energy has saturated at a level of 
about 1 per cent of the initial kinetic energy. Both fields 
oscillate irregularly after that, and are weakly turbulent. 
It is unclear how such a small magnetic excitation suc- 
ceeds at shutting down such a large fraction of the flow. 
As will be shown later, this large drop is associated with 
the instability of the large scale flow. The inset shows 
the full time history of Ey and Em, for Ry — 3276.6 
and Rm = 1092.2 when the turbulence is fully devel- 
oped. The dashed line illustrates, for comparison, how 
simply the magnetic energy exponentiates and saturates 
in the laminar steady-state regime (Ry = 62.7). Figure 
shows in detail the suppression of the flow, manifested 
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FIG. 7: Time history of the total energy in the dynamo sim- 
ulation with Ry = Rm ~ 416.2. The shaded area is a blow 
up of the shaded region in Fig. and corresponds to the 
hydrodynamic oscillations. 



as a drop in the total energy, at t m 70. 

These oscillations between the hydrodynamic laminar 
and turbulent regime in the Roberts flow have been pre- 
viously found by Feudel et al. ||. The authors pointed 
out that in this regime, close to the threshold R C M the 
dynamo exhibits an intermittent behavior, with bursts 
of activity. The oscillatory flow is stable to small per- 
turbations (e.g. numerical noise in the code), but as 
the magnetic energy grows the flow is perturbed by the 
Lorentz force and goes to a weakly turbulent regime. As 



noted in Ref. [8|, if Rm is close to R C M then the mag- 
netic field decays, the flow relaminarizes and the process 
is repeated. However, as observed in Fig. [f)| if Rm is 
large enough the weakly turbulent flow can still excite a 
dynamo, and the magnetic field keeps growing exponen- 
tially until reaching the non-linear saturation even after 
the hydrodynamic instability takes place. 

Figure [5fa) shows the temporal growth of several 
Fourier components of the magnetic field in the laminar 
regime (Ry = 62.7). A straightforward exponentiation, 
followed by a flat, steady-state, lcvcling-off exhibits the 
same growth rate for all harmonics. This indicates the ex- 
istence of a simple unstable normal mode which saturates 
abruptly near t ~ 180. The behavior is much noisier for 
R v = 416.2 and 3276.6 as shown in Figs. |SJ» and^c). 
Note that in the simulation with Ry = 416.2, for t < 70 
all the magnetic modes oscillate with the same frequency 
as the hydrodynamic oscillations. In Figure [HI the dot- 
ted line and solid line above are, respectively, for k = 1 
and k = 2. The remaining four are for k = 9 through 
11. The modes in between occupy the open space in be- 
tween more or less in order. The same modes are shown 
for Ry = 416.2 in Fig. OJc), which illustrates a broad 
sharing of Eb among many modes and a consequent ex- 
citation of small-scale magnetic components. 

Plots of the kinetic and magnetic fields are shown in 
Figs. [UJ The left column shows the velocity field in the 
saturated state for Ry = 62.7, and the right column 
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FIG. 8: Evolution of the magnetic energy in different shells in 
Fourier space: (a) Rv — 62.7 and Rm = 78.4 (laminar flow), 
(b) Rv = Rm = 416.2 (periodic case), (c) Rv = 3276.6 and 
Rm — 1092.2 (turbulent regime) . The dotted line corresponds 
to k = 1, solid line to k = 2, and the dashed lines to k = 
9,10,11,12. 



shows the magnetic field at the same time. The arrows in- 
dicate the vector components in the planes shown and the 
colors indicate the strengths of the perpendicular compo- 
nents. Figures E^a) and (b) are for the plane z = and 
Figs. EJc) and (d) are for the plane y = 0. Figs. OJe) and 
(f) are for the plane y = tt/2. The velocity configuration 
shown in Fig. Ela) is quite similar to the way it looks at 
t = 0, but the z-dependences apparent in Figs. He), (d), 
and (f) are not present in the initial flow. 

Figs. EJare similar color plots for the saturated regime 
for Ry = 416.2. All the same quantities are being dis- 
played at the same planes as in Figs. El The initial con- 
ditions are no longer recognizable in the saturated state, 
but is not yet sufficiently disordered that one would be 
forced to call it "turbulent" . Moreover, note that the four 
"cells" characteristic of the laminar Roberts flow [Fig. 
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FIG. 9: (Color online) Plots of the kinetic and magnetic fields for the saturated regime of the run with Rv = 62.7 and 
Rm = 78.4: (a) cut at z = 0, vz in color and vx, vy indicated by arrows, (b) same as in (a) for the magnetic field, (c) cut at 
y — 0, vy in color and vx, vz indicated by arrows, (d) same as in (c) for the magnetic field, (e) same as in (b) but for a cut at 
y — 7r/4, and (f) same as in (e) for the magnetic field. 



EUa)] are not present in this late stage of the dynamo. 
During the early kinematic regime, when the hydrody- 
namic oscillations are observed, a slightly deformed ver- 
sion of these cells can be easily identified in the flow (not 
shown). When the magnetic energy grows due to dynamo 
action, the flow is unable to maintain this flow against 
the perturbation of the Lorentz force. This causes the 
large-scale flow to destabilize, and the kinetic energy in 
the shell k = 1 drops by a factor of two. This instability 
of the large-scale modes is associated with the large drop 
of the kinetic and the total energy at t ss 70 (Fig. [7J. 

By contrast, the same fields are exhibited in the same 
planes in Figs. ^2 m the saturated regime for Ry = 
3276.6. Here the truly turbulent nature of the flow is now 
apparent, particularly in the highly disordered magnetic 
field plots in the right-hand column. 



Figurelr21is a three-dimensional perspective plot of the 
kinetic and magnetic energy density for Ry = 62.7 at a 
late time in the saturated regime. The kinetic energy 
distribution (on the left) is not much different than it 
was at t = 0. The helical properties of the Roberts flow 
can be directly appreciated in the field lines indicated in 
black. In this regime, the flow is still laminar as pre- 
viously indicated. The magnetic field is stretched and 
magnetic energy amplified in the four helical tubes, and 
then expelled out of the vortex tubes, accumulating in 
the stagnation points 0, Q . Since the velocity field has 
no dependence in the z-direction, the magnetic field that 
can be sustained by dynamo action has to break this sym- 
metry and displays a clear periodicity in this direction. 
The same energy densities are exhibited at a late time 
for the case of Ry = 3276.6 in Fig. and the highly 
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FIG. 10: (Color online) Plots of the kinetic and magnetic fields 
Labels and fields are as in Fig. |3] 

filamented and disordered distributions characteristic of 
the turbulent regime are again apparent. Note however 
that still some helicity can be identified in the velocity 
field lines shown. 

In Ref. a suppression of small scale turbulent fluc- 
tuations and an evolution of the system to a state with 
effective magnetic Prantdl number of order one was ob- 
served in the nonlinear saturation of the turbulent dy- 
namo. Here a similar effect is observed, although the 
suppression of small scales is weaker probably due to the 
presence of the external forcing at k ps 1 which does not 
leave room for a large scale magnetic field to develop. 
Figure ^] shows the time evolution of the kinetic and 
magnetic energy spectra in the run with Ry = 3276.6 
and Rm = 1092.2. While at early times the magnetic 
energy spectrum peaks at small scales (k « 9), at late 
times the magnetic spectrum is flat for small k and drops 
together with the kinetic energy. The kinetic spectrum 
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for the saturated regime of the run with Rv = Rm = 416.2. 



is strongly quenched and has a large drop at small scales. 



IV. SUMMARY AND DISCUSSION 

One apparent outcome of these computations has been 
to confirm the intuitive impression that dynamo amplifi- 
cation of very small magnetic fields in conducting fluids 
is easier if mechanical helicity is present. This is true in 
velocity fields which arc both turbulent and laminar. The 
values of R C M which are the lowest found (~ 10) are well 
below those in several existing experimental searches. 

It is also somewhat reassuring to find that the qual- 
itative behavior of dynamo thresholds with decreasing 
viscosity (increasing Reynolds number at fixed U) is as 
similar as it is to that found for the non-helical TG 
flow in Ref. 0. In particular, since the simulations 
discussed here were forced at almost the largest scale 
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FIG. 11: (Color online) Plots of the kinetic and magnetic fields for the saturated regime of the run with Rv = 3276.6 and 
Rm = 1092.2. Labels and fields are as in Fig. |U] 




FIG. 12: (Color online) Visualization of the kinetic (left) and FIG - 13: ( Color online) Visualization of the kinetic (left) and 

magnetic energy density (right) for the saturated regime of magnetic energy density (right) for the saturated regime of 

the run with R v = 62.7 and R M = 78.4. Velocity field lines the run with R V = 3276.6 and R M = 1092.2. Velocity field 

are indicated in black. lines are indicated in black. 



was reached usi ng only DNS, while for the TG flow two 
available in the periodic domain, a turbulent regime for different models [13,13 f° r the small scales were needed. 
Pm < 1 where R C M is approximately independent of Pm The similarities in the behavior of the threshold for the 
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That being said, it should be admitted that the Roberts 
flow in a way exhibits a richer set of possibilities in that 
the dynamo activity is somewhat different in each of the 
three regimes (laminar and steady-state, oscillatory, and 
turbulent). 

Dynamo action is to be regarded as of many types 
and situation-dependent. The forms of the magnetic 
fields developed and their characteristic dimensions are 
determined to a considerable extent by the mechanical 
activity that excites them and by the geometric setting 
in which they take place. If it is desired to apply the 
theoretical and computational results to planetary dy- 
namos or laboratory experiments, then rectangular peri- 
odic conditions appear to be a constraint that should be 
dispensed with as soon as feasible. 



FIG. 14: (Color online) Kinetic [thick (blue) line] and mag- 
netic energy spectra for different times for the simulation with 
Rv = 3276.6 and R M = 1092.2. 

two flows for Pm small enough brings more confidence 
to the ability of subgrid scale models of MHD turbulence 
to predict results in regimes of interest for astrophysics 
and geophysics that are today out of reach using DNS. 
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